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GLOSSA.Y  OF  TERMS 


All  symbols,  subscripts,  and  nondimensionalized  hydrodynamic  coef- 
ficient expressions  appearing  in  the  body  of  this  report  are  defined 
below.  Any  dimensions  involved  will  be  consistent  with  the  foot-pound- 
second  system  of  units.  All  angles  are  in  degrees. 


SYMBOLS 

Symbol 

AR 

B 


b 

c 

c 

c 

r 


1 

a 


I 

x 


T 

‘z 


K 


k 


1 


Definition 

2 

Aspect  ratio  of  surface  = b /s. 

Buoyancy  force  which  is  positive  upwards,  B * pV  (except 
where  used  in  fluid  kinetic  energy  equations) . 

Exposed  root- to- tip  span  of  surface. 

Chord. 

Mean  chord  of  surface  = s/b  - b/AR. 

Root  chord  of  surface. 

Tip  chord  of  surface. 

Added  moment  of  inertia. 

Moment  of  inertia  of  an  underwater  vehicle  about  the 
x-axis. 

Added  moment  of  inertia  about  the  x-axis. 


Moment  of  inertia  of  an  underwater  vehicle  about  the 
y-axis. 

Moment  of  inertia  about  the  y-axis  of  the  fluid 
displaced  by  the  body. 

Moment  of  inertia  of  an  underwater  vehicle  about  the 
z-axis. 

Moment  about  the  x-axis  exerted  on  the  body  by  the 
dynamic  pressure  of  the  surrounding  fluid. 

Lamb's  inertial  coefficient  for  a prolate  ellipsoid 
in  axial  flow. 
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Symbol  (Continued)  Definition 

1<2  Lamb's  inertial  coefficient  for  a prolate  ellipsoid 

in  cross  flow. 


k'  Coefficient  for  added  moment  of  inertia  for  a flat  plate. 

P 

k^  Lamb's  coefficient  for  added  moment  of  inertia  for  a 

prolate  ellipsoid  body. 

S,  Overall  length  of  the  vehicle. 

M Moment  about  the  y-axis  exerted  on  the  body  by  the 

dynamic  pressure  of  the  surrounding  fluid. 

m Mass  of  the  vehicle  including  the  water  in  the  free- 

flooding  spaces. 

m Additional  mass. 

a 

m,,.  Mass  of  fluid  displaced  by  hull  envelope. 

dr 

N Moment  about  the  z-axis  exerted  on  the  body  by  the 

dynamic  pressure  of  the  surrounding  fluid. 

p Component  of  angular  velocity  about  th  i body 

fixed  x-axis. 

q Component  of  angular  velocity  about  the  body 

fixed  y-axis. 

r Component  of  angular  velocity  about  the  body 

fixed  z-axis. 

s Exposed  planform  area  of  surface. 

s,  Maximum  cross-sectional  area  of  the  body. 

b 

T Total  kinetic  energy  of  the  body-fluid  system. 

Tf  Kinetic  energy  of  the  fluid  representing  the  system  cZ 

impulsive  pressures  exerted  by  the  surface  of  the  bodj’ 
on  the  fluid  during  acceleration. 

Tveh  Kinetic  energy  of  the  vehicle  representing  linear  and 
4'!  angular  momentum  of  the  body. 

t Time. 

U Linear  velocity  of  origin  of  body  axes  relative  to  a 

fluid- fixed  axis  system. 

ii 
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Symbol  (Continued) 


Definition 


u 

V 

v 

w 

X 


neb 


neg 


nht 


cb 


bp 


Component  of  U along  the  body  x-axis. 

Displaced  volume  of  the  hull  envelope. 

Component  of  U along  the  body  y-axis. 

Component  of  U along  the  body  z-axis. 

Force  along  the  x-axis  exerted  on  the  body  by  the 
dynamic  pressure  of  the  surrounding  fluid. 

x-distance  from  nose  to  cb. 

x-distance  from  nose  to  eg. 

x-distance  from  nose  to  centroid  of  the  horizontal 
tail  fins. 

Longitudinal  axis  of  the  body  fixed  coordinate  axis 
system. 

Force  along  the  y-axis  exerted  on  the  body  by  the 
dynamic  pressure  of  the  surrounding  fluid. 

Transverse  axis  of  the  body  fixed  coordinate  axis 
system. 

Force  along  the  z-axis  exerted  on  the  body  by  the 
dynamic  pressure  of  the  surrounding  fluid. 

Vertical  axis  of  the  body  fixed  coordinate  axis 
system. 

The  z-distance  to  the  center  of  buoyancy  (cb)  from  the 
center  of  gravity  (eg) , positive  down. 

Angle  of  attack. 

Angle  of  side-slip. 

Deflection  angle  of  bowplane  (or  sailplane) . 

Deflection  angle  of  rudder. 

Deflection  angle  of  sternplane. 
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Symbol  (Continued) 


Definition 


X 

w 

* 

♦ 

P 

0 

J 


Planform  taper  ratio  of  surface  ■ ct/cr* 

Angular  velocity. 

Yaw  angle. 

Roll  angle. 

Mass  density  of  seawater. 

Pitch  angle. 

The  component  in  the  plane  of  a surface  of  the  perpen- 
dicular distance  between  the  axis  of  rotation  and  the 
centroid  of  the  area  of  the  surface. 

A prime  over  any  symbol  signifies  nondimensionalization. 

A dot  over  any  symbol  signifies  differentiation  with 
respect  to  time. 


SUBSCRIPTS 


Subscript  Definition 

a Added  or  additional  inertia., 

b Body, 

bp  Bow  plane, 

cb  Distance  from  eg  to  cb. 

eg  Distance  to  eg. 

df  Displaced  fluid, 

f Fluid, 

ht  Horizontal  tail. 

Jlvt  Lower  vertical  tail  section. 
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Subscript  (Continued)  Definition 

n(  ) Referenced  from  nose  to  some  point,  ( ),  e.g 
implies  nose-to-cb. 

p Plate. 

r Rudder . 

s Sternplane. 

sh  Shroud. 

t Tip. 

uvt  Upper  vertical  tail  section, 

veh  Vehicle. 

x Pertaining  to  the  x-axis. 

y Pertaining  to  the  y-axia. 

z Pertaining  to  the  z-axis. 
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NONDIMENS TONAL I ZED  HYDRODYNAMIC  COEFFICIENT  EXPRESSIONS 


V 1 

Xu 

yt  m . 

X. 

V 

X . 
u 

" 75)7F  ’ 

V 

TWF 

v 1 

V*  - 

Y. 

V 

Y . 
u 

“ WpP  » 

I • ■ 

V 

r*  1 

Zu 

7 1 ■ 

z. 

V 

Z . 
u 

■ (>l)pt»  • 

4 • * 
V 

<*s)pRa 

vt 

. 

y t - 

. *A 

x* 

x 4 

Ci)pV 

V f 

YJ>  . 

V 1 . 

.ya 

Y’  . 
p 

I . ■ 

q 

(4)  Pi4 

ZP  . 

7 1 m 

z. 
q , 

z . 

p 

Os)pF 

4 • * 

q 

Oj)pF 

K'u 

K. 

u 

K'.  - 

V 

K. 

V 

* '(*s)pF  * 

<4)pF 

M'  . 
u 

M. 

u 

M’ . - 

V 

M. 

V 

" O^F  ’ 

(%)p'F 

N’  . 
u 

N. 

u 

N\  - 

V 

N. 

V 

■ • 

Os)pF 

v V 

Kf, 

K\  - 

q 

*4 

K . 
P 

WpF  * 

Os)pF 

MV 

4MU 

T5)pF  * 

(4>p*° 

N'  ^ 

P 

N. 

> P 

s'$  ■ 

N. 

a_ 

'VppF’ 

15)pF 

Note  that  X’u  ■ 9 X'/Su’,  X'p  ■ 9X'/9p'»  etc. 
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INTRODUCTION 


BACKGROUND 

This  report,  is  primarily  a compilation  of  methods  for  analytically 
predicting  the  acceleration  hydrodynamic  coefficients  from  the  geometric 
and  mass  distribution  characteristics  of  an  underwater  vehicle.  It  rep- 
resents the  current  results  of  an  effort  to  fin  ! the  most  generally  ap- 
plicable and  accurate  techniques  for  calculating  these  coefficients. 

The  need  for  such  analytical  techniques  grew  from  the  many  unique 
underwater  vehicle  requirements  inherent  to  the  Naval  Coastal  Systems 
Laboratory  mission  areas.  Typically  these  requirements  include  design, 
modification,  or  simulation  of  towed,  “ethered,  or  free  vehicles  for 
which  analytical  calculation  of  the  hv  .odynamic  coefficients  is  neces- 
sary. This  necessity  arises  from  the  need  for  some  combination  of  the 
following: 

1.  Flexibility.  The  user's  project  is  in  a design  phase 
where  vehicle  requirements  imposed  by  hardware  sub- 
systems are  subject  to  change  and  it  is  desirable  to 
assess  the  effect  of  many  more  vehicle  configurations 
than  would  be  practical  for  experimental  determination. 

2.  Responsiveness.  The  project  is  of  a research  and  develop- 
ment nature  such  that  there  is  a requirement  for  rapid  and 
frequent  updating  of  the  vehicle's  predicted  performance 
in  response  to  changes  in  performance  goals  and/or  subsys- 
tem hardware. 

3.  Economy.  Cost  constraints  on  the  project  will  not  permit 
use  (or  extensive  use)  of  model  basin  testing,  but  the 
project  still  requires  pome  degree  of  assurance  that  there 
are  no  potential  catastrophic  failures  due  to  problems 
such  as  dynamic  instability. 

To  more  effectively  meet  each  of  these  requirements,  all  of  the 
methods  presented  herein,  as  well  as  methods  for  computing  all  of  the 
remaining  hydrodynamic  coefficients,  have  been  computer- Implemented  in 
the  complete  stability  and  control  analysis  program  GEORGE.  The  GEORGE 
computer  program  uses  vehicle  geometric  and  mass  distribution  information 
to  compute  vehicle  dynamic  characteristics. 
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It  should  be  pointed  out  that  the  procedure  yielding  the  highest 
degree  of  confidence  for  an  underwater  vehicle  design  involves  the  use 
of  both  model  basin  testing  and  analytical  techniques.  An  example  of 
a desirable  sequence  of  design  steps  would  be  the  Implementation  of  a 
computerised  hydrodynamic  math  model  of  the  proposed  vehicle  using  con- 
ceptual information  available  in  the  very  early  design  phases.  Using 
this  model  and  incorporating  updated  information  as  it  becomes  avail- 
able, the  effect  of  various  sizing,  distributing,  and  shaping  ideas  on 
the  vehicle  dynamics  can  be  assessed. 

As  soon  as  a reasonably  firm  design  configuration  is  laid  out, 
model  basin  tests  can  be  performed  to  validate  the  computer  modeling. 
These  tests  can  be  much  less  extensive  than  those  necessary  to  produce 
the  parametric  information  required  as  a basis  for  new  design  tech- 
niques— techniques  sufficiently  general  to  maintain  predictive  accuracy 
for  the  redesigns  or  modifications  frequently  necessary  in  vehicle  de- 
sign. This  reduction  in  model  testing  requirements  is  primarily  the 
result  of  having  a computer  model  to  use  in  defining  what  needs  to  be 
verified  experimentally,  and  in  providing  a sophisticated  tool  for 
describing  complicated  trends,  thus  reducing  data  requirements. 

The  validated  computer  simulation  model  can  then  be  used  in  many 
ways  to  support  the  project.  Trajectory  simulations  can  be  made  for 
given  operational  conditions  so  that  motion  information  can  be  genera- 
ted for  performance  evaluation  or  definition  of  hardware  specifications. 
Examples  of  work  performed  using  simulations  in  this  manner  are  trajec- 
tory analyses  of  submarine- launched  acoustic  countermeasure  devices  and 
motion  prediction  of  a swimmer  delivery  vehicle  to  assess  the  motion  in- 
fluence on  the  design  of  a high-resolution,  side- looking  sonar.* 

Another  use  has  been  to  provide  the  hydrodynamic  coefficients  for 
analog  or  hybrid  simulation  models  to  study  man-machine  interface  prob- 
lems or  to  provide  training.  An  example  of  recent  work  in  this  area  is 
the  implementation  of  a hybrid  model  of  the  Mine  Neutralization  Vehicle 
(MNV)  as  a performance  evaluation  and  training  device.  The  MNV  simula- 
tion solves  the  nonlinear  equations  of  motion  for  a tethered  vehicle 
(including  four  propulsors,  bowplanes,  and  cable  tension)  and  produces 
a motion  display  for  either  body  or  inertial  reference  frame,  all  in 
real  time.  Of  course,  the  computer  simulation  model  is  most  often  used 
to  assess  the  influence  of  design  variations  of  the  vehicle  hydrodynamic 
and  control  system  parameters  on  dynamic  stability  and  control 
characteristics. 


*Trajectory  simulations  have  usually  been  performed  using  suitable  modi- 
fied versions  of  the  standard  six-degrees-of-freedom  submarine  simulation 
computer  program,  ZZMN,  developed  by  DTNSRDC. 
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Having  briefly  stated  the  background  for  the  prediction  and  Simula' 
tion  capability  regarding  the  stability  and  control  of  underwater 
vehicles  at  NCSL,  attention  will  now  be  focused  on  the  analytical  tech- 
niques for  assigning  a value  to  acceleration,  nonlinear,  and  cross- 
coupled  hydrodynamic  coefficients.*  These  techniques  have  been  applied 
to  a wide  variety  of  body  types,  including  torpedoes,  submarines,  swim- 
mer delivery  vehicles  (SDV's),  and  many  other  special  purpose  vehicles 
such  as  the  Deep  Submergence  Rescue  Vehicle  (DSRV) , the  Shipboard  Mine- 
hunting System  (SMS)  towed  sonar  body,  and  submarine- launched  torpedo 
countermeasures  devices. 

Generality  in  applicability  has  been  sought;  that  la,  use  of  a 
method  that  yields  accurate  results  for  only  a specific  class  of 
vehicles,  such  as  torpedoes,  has  been  avoided.  A more  generally  appli- 
cable and,  in  the  experience  of  the  authors,  more  accurate,  formulation 
based  on  theoretical  derivation  and  experimental  verification  is  pre- 
sented. The  relative  importance  of  each  hydrodynamic  coefficient  as  it 
contributes  to  the  responsive  motion  of  the  vehicle  is  noted.  In  addi- 
tion, important  relationships  between  acceleration,  cross-coupled,  and 
nonlinear  hydrodynamic  coefficients  are  derived  from  potential  flow 
theory. 


ACCELERATION  COEFFICIENTS 

The  forces  represented  by  acceleration  hydrodynamic  coefficients 
are  attributed  to  the  acceleration  of  fluid  particles  surrounding  a sub- 
merged body.  Some  of  the  forces  appear  in  the  equations  of  motion  as 
additions  to  terms  representing  the  inertia  of  the  body.  They  are  known 
as  aided  mass  or  msss  accession  terms  and,  when  combined  with  the  in- 
ertial parameters  of  a body,  have  been  called  apparent,  virtual,  and 
hydrodynamic  mass. 

Analytical  methods  for  estimating  the  acceleration  coefficients  of 
the  bare  hull  differ  from  the  methods  used  for  the  contributions  of  ap- 
pendages in  that  they  are  not  estimated  using  semi-empirical  formulae. 
Rather,  the  body  is  assumed  to  be  moving  through  an  infinite,  inviscid, 
circulation-free  fixed  fluid.  For  these  conditions,  Lamb^x)  derives  the 
complete  expression  for  the  kinetic  energy  of  the  body-fluid  system 
which,  coupled  with  inertia  coefficients  for  prolate  ellipsoids^ , 
yields  estimates  of  body  bare-hull  acceleration  derivatives. 


*Although  techniques  are  also  available  for  the  prediction  of  static 
and  dynamic  velocity  terms  with  viscous  effects  included,  the  scope 
of  this  report  is  limited  to  the  acceleration  coefficients. 

(1^Lamb,  Sir  Horace,  Hydrodynamics,  Dover  Publications,  New  York, 

Sixth  Edition,  1932.  pps  172  and  154. 
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The  basis  for  the  added  masses  and  moments  due  to  body  appendages, 
such  as  fins  and  bovplanes,  is  largely  semi-empirical.  Theoretical  ex- 
pressions for  thin,  flat  platea  moving  in  a perfect  fluid  have  been 
derived  from  hydrodynamic  theory  with  experimentally  determined  correc- 
tion factors  for  non-ideal  conditions  and  shapes.  A great  deal  of  the 
primary  work  done  in  this  area  was  accomplished  by  Munk(a)(3),  with 
later  significant  experimental  work  by  Gracey^,  and  Malvestuto  and 
Gale(B). 


ACCELERATION  COEFFICIENTS  IN  POTENTIAL  FLOW  THEORY 


The  forces  and  moments  represented  by  the  acceleration  hydrodynamic 
coefficients  can,  to  a very  great  extent,  be  modeled  as  potential  flow 
phenomena.  Neglecting  the  details  of  the  boundary  layer  In  modeling 
acceleration-dependent  forces  and  moments  acting  on  a submerged  body 
yields  quite  satisfactory  results  for  most  stability  and  control  simula- 
tion. Coefficient  contributions  due  to  appendages  are  most  noticeably 
influenced  by  boundary  layer  effects.  Acceleration  coefficient  values 
predicted  using  potential  flow  theory  are  probably  most  suspect  in  their 
application  to  fins  that  have  a significant  portion  of  their  planform 
area  in  separated  flow.  There  are  a few  considerations  concerning  these 
sources  of  error  that  should  be  mentioned  before  launching  into  further 
development . 

Acceleration-dependent  forces  can  arise  from  changes  in  the  bound- 
ary layer,  resulting,  in  turn,  from  changes  in  angle-of-attack  or  drift 
angle.  There  are  several  interesting  facts  in  this  regard.  One  is  that 
the  speed/power  relationship  for  underwater  vehicles  generally  encourages 
a designer  to  fair  his  vehicle  so  that  the  boundary  layer  has  a signifi- 
cant thickness  over  as  little  of  the  vehicle  as  possible.  The  result 


(^National  Advisory  Committee  for  Aeronautics  T.N.  No.  197,  Sane  Tables 
of  the  Factor  of  Apparent  Additional  Mass,  by  Max  M.  Munk,  1924. 

(3)Munk,  Max  M. , Fluid  Mechanics,  Part  II,  Aerodynamic  Theory,  Vol.  1, 
edited  by  W.  F.  Durand,  Dover  Publications,  Inc. , New  York,  1934. 

National  Advisory  Committee  for  Aeronautics  T.N.  No.  707,  The  Addi- 
tiona 1-Mass  Effect  of  Plates  as  Determined  by  Experiments , by 
William  Gracey,  1941. 

(^National  Advisory  Cousnittee  for  Aeronautics  T.N.  No,  1187,  Formulas 
for  Additional-Mass  Corrections  to  the  Moments  of  Inertia  of 
Airplanes,  by  Frank  S.  Malvestuto,  Jr.  and  Lawrence  J.  Gale,  1946. 
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is  that  the  pressure  distribution  over  80  percent  of  the  surface  area 
of  a typical  underwater  vehicle  can  be  described  very  accurately  by 
potential  flow  theory.  Consequently,  the  acceleration  coefficients 
most  important  to  the  description  of  underwater  vehicle  dynamics  can 
usually  be  accurately  computed  as  potential  flow  phenomena.  Some  less 
important  acceleration  coefficients,  especially  those  resulting  from 
the  imbalance  of  added  mass  effects  distributed  fore  and  aft  of  the  eg, 
can  be  significantly  erroneous  if  the  boundary  layer  is  neglected.  Un- 
fortunately, those  coefficients  influenced  by  boundary  layer  effects 
suffer  from  similar  problems  in  experimental  model  testing;  e.g.,  im- 
proper Reynolds  number  scaling  leads  to  nondynamic  similarity  in  the 
modeled  boundary  layer. 

Describing  acceleration-dependent  effects  on  fins  is  simplified 
somewhat  by  the  vehicle  designer's  avoidance  of  the  sizing  and  place- 
ment of  fins  so  that  they  lose  effectiveness  in  regions  of  separated 
flow.  Consequently,  more  often  than  not,  treatment  of  the  added  mass 
effects  of  fins  as  if  they  are  in  unspoiled  flow  provides  good  results. 
Tube- launched  devices  are  often  exceptions  to  this  rule.  They  are  gen- 
erally constrained  by  tail  size  and  after-body  taper  so  that  large  por- 
tions of  the  fins  are  in  separated  flow. 

Decisions  on  "effective"  fin  aspect  ratios  based  on  past  experi- 
ence with  similar  vehicles  may  be  necessary  to  achieve  ai  curacy  in 
those  coefficients  with  significant  fin  contributions.  Since  the  fins 
on  a torpedo  are  generally  small,  their  contributions  to  the  major  coef- 
ficients of  importance  are  less  significant  than  is  the  case  with  other 
types  of  underwater  vehicles.  Roll  coefficients  are  important  ex- 

ception since  they  are  due  almost  completely  to  fin  cu?vi i ibutions . 

The  acceleration,  cross-*coupled,  and  nonlinear  coefficient  equali- 
ties derived  in  Appendix  A must  be  used  with  great  caution  when  simu- 
lating vehicle  trajectories.  The  NCSL  approach  has  been  to  assess  the 
importance  of  coefficient  accuracy  using  two  computer- implemented  tools. 
One  computer  program  performs  root  locus  sensitivity  studies  for  the 
linear  coefficients  on  the  vehicle  of  interest.  The  required  accuracy 
is  determined  on  a coefficient-by-coefficient  basis,  taking  advantage 
of  the  insight  that  can  be  gained  by  a control  system  analysis  approach 
to  the  vehicle  system. 

A more  thorough  evaluation  of  the  math  model  simulating  the  vehicle 
fluid  system,  including  the  nonlinear  and  cross-coupled  terms,  must  be 
done  on  a maneuver-by-maneuver  basis.  To  assess  the  relative  importance 
of  each  contributing  force  and  moment,  the  authors  have  modified  the 
standard  six-degrees-of-freedom  submarine  simulation  computer  program^ 


^6)Naval  Ship  Research  and  Development  Center  Test  and  Evaluation  Report 
P-433-H-01,  User's  Guide  NSRDC  Digital  Program  for  Simulating  Subma- 
rine Motion  ZZMN-Revision  1.0,  by  Ronald  W.  Richards,  June  1971. 
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to  plot  the  following  on  any  maneuver:  (1)  time  histories  of  the  force 
(or  moment)  of  each  contributing  coefficient,  propulsor,  buoyancy/ 
gravity,  and  control  surface  term  together  with  the  total  force  (or 
moment);  and  (2)  time  histories  of  the  motions  of  the  vehicle  with  the 
contribution  of  a single  term  removed,  and  repeated  for  each  term.  Know- 
ing the  degree  of  accuracy  expected  in  the  computed  value  for  each  c ef- 
ficient, one  can  analyze  the  plots  on  the  maneuvers  of  interest  and  assign 
a level  of  confidence.  This  procedure  will  give  an  indication  of  which 
coefficients  are  most  critical  in  achieving  accurate  simulation  on  a 
given  maneuver. 

Using  these  techniques  and  a multitude  of  comparisons  with  planar 
motion  mechanism  and  rotating  arm  results,  the  methods  presented  in  this 
report  are,  on  the  basis  of  the  authors'  experience,  considered  to  be 
quite  accurate  on  any  hydrodynamically-f aired  vehicle  for  all  but  ex- 
treme maneuvers.  This  is  especially  true  for  linear  models,  as  will  be 
shown  by  examples  in  Appendix  B. 

The  modeling  accuracy  achieved  for  acceleration  terms  under  these 
assumptions  is  quite  contrary  to  the  extremely  poor  results  obtained 
when  the  same  assumptions  are  applied  to  velocity-dependent  terms.  The 
assumptions  are:  (1)  the  body  is  submerged  so  that  there  are  no  near- 
surface or  near-bottom  effects;  (2)  the  fluid  is  inviscid;  and  (3)  there 
is  zero  circulation. 


1 


i 


I 


i 

\ 
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Under  these  conditions  the  momentum  equations  for  the  body-fluid 
system  may  be  written  in  the  six  Lagrange  equations: 
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the  body  by  the  pressure  of  the  surrounding  fluid,  and  the  total  kinetic 
energy  T may  be  written 
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with  Tygjj  representing  the  linear  and  angular  momentum  of  the  vehicle 
itself,  and  Tf  representing  the  system  of  impulsive  pressures  exerted 
py  the  surface  of  the  solid  on  the  fluid  in  the  supposed  instantaneous 
generation  of  motion  from  rest. 

If  we  rewrite  Equations  (1)  and  (2),  isolating  the  terms  due  to  Tf, 
we  obtain  expressions  for  the  forces  exerted  on  the  body  by  the  dynamic 
pressure  of  the  surrounding  fluid,  e.g.: 


X - 

and 


K - 

etc. 

The  fluid  kinetic  energy,  Tf,  can  be  expressed  as  a quadratic  form 
of  the  body  axis  velocities  u,  v,  w,  p,  q,  and  r with  no  explicit  time 
dependence.*  The  quadratic  expression  for  Tf  may  be  written^: 

2Tf  * Au2  + Bv2  + Cw2  + 2A’vw  + 2B'wu  + 2C'vu  + Pp2  + Qq2 

+ Rr2  + 2P'qr  + 2Q'rp  + 2R'pq  + 2Lup  + 2Mvq  + 2Nwr 

+ 2F(vr  + wq)  + 2G(wp  + ur)  + 2H(uq  + vp)  + 2F'(vr  - wq) 

+ 2G' (wp  - ur)  + 2H'(uq  - vp) . (4) 

The  body  axes  are  generally  chosen  so  that  the  xz-plane  is  a plane 
of  symmetry  (port-starboard  symmetry) . This  eliminates  nine  of  the 
coefficients.  A*,  Cf,  P',  R’ , L,  M,  N,  G,  and  G’,  since  by  symmetry  argu- 
ments, sign  changes  in  v,  p,  or  r should  not  change  Tf.  Letting 

F + F'  - Flt  F - F*  - F2,  H + H*  - Hx,  and  H - H’  - H2,  (5) 


(1)ibid. 

*Recent  investigations  indicate  that  there  can  be  significant  time  vari- 
ant terms  arising  from  vortices  shed  from  the  forward  appendages  and 
hull.  The  mathematical  modeling  of  this  phenomenon  is  accomplished 
using  force  and  moment  terms  in  the  equations  of  motion  that  are  convo- 
luted in  time.  These  terms  are  outside  the  set  of  acceleration,  non- 
linear, and  cross-coupled  hydrodynamic  coefficients  under  consideration 
here. 


!Ll_ 

dt  3u  r 3v  ^ 


3Tf 

3w 
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dt  3p  w 3v 


3Tf  9Tf  JT 
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the  remaining  12  coefficients  are  included  in  the  equation 
2Tf  “ Au2  + Bv2  + Cw2  + 2B'wu  + Pp2  + Qq2  + Rr2 

+ 2Q'rp  + 2F^vr  + 2F2wq  + 2H^uq  + 2H2vp  . (6) 

The  complete  expressions,  assuming  no  planes  of  symmetry,  are  pre- 
sented in  Appendix  A. 


Proceeding  as  indicated  in  Equation  (3),  the  full  six-degrees-of- 
freedom  momentum  equations  can  be  written  in  terms  of  the  12  remaining 
coefficients*7*.  The  coefficients  are  written  in  hydrodynamic  coeffi- 
cient notation  as  follows. 


Added  inertia  terms: 

A = -X.,  P = -K.,  B = -Y.,  (7) 

u p’  v* 

Q - -M.,  C - -Z.,  and  R = -N.  . (7) 

q w r 

Inertial  intermodal  coupling  terms: 

FX  = ~Yi  " “V  F2  - -Z.  = -M*,  and  H2  - -Y.  = -K^  . (8) 

The  remaining  three  are 

Bf  = -X^  = -Z., , Q»  = -N.  = -Kf,  and  ^ = -X.  = -M&  . (9) 

If,  as  shown  in  the  complete  expressions  in  Appendix  A,  there  is 
also  fore-aft  symmetry  (symmetry  about  the  yz-plane) , then,  in  addition 
to  the  aforementioned  zeroed  terms 

B'  = Q'  = Fx  = F2  = 0 . 

Vehicles  having  symmetry  about  both  the  xz-plane  and  the  xy-plane  have 
the  additional  zeroed  terms 

B'  = Q*  = = H2  = 0 . 


Massachusetts  Institute  of  Technology,  Instrumentation  Laboratory 
Report  R-570-A,  Deep  Submergence  Rescue  Vehicle  Simulation  and  Ship 
Control  Analysis,  by  Charles  Broxmeyer,  Pierre  P.  Dogan,  et  al, 
1967,  p.  31. 
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These  relations  will  be  used  extensively  in  the  following  coeffi- 
cient formulations.  It  should  be  noted  that  the  coefficient  identities, 
such  as  Equations  (8)  and  (9),  result  from  derivation  under  no  assump- 
tions of  symmetry.  The  existence  of  symmetry  results  in  some  of  •‘he 
terms  being  equal  to  zero.  The  complete  set  of  identities  of  accelera- 
tion, nonlinear  and  cross-coupled  coefficients  is  presented  in  Appendix  A. 


METHODS  FOR  COMPUTING  ACCELERATION  HYDRODYNAMIC  COEFFICIEN 


INTRODUCTION 

Methods  for  estimating  the  acceleration-dependent  hydrodynamic  co- 
efficients are  given  in  this  section.  Each  coefficient  is  described 
using  the  following  formats 

1.  Introductory  description. 

2.  Relative  importance  - major  or  minor,  and  under  what 
circumstances . 

3.  Causal  relationship  - important  assumptions. 

4.  Analytical  expression  - including  normal  range  of  coefficients 
in  the  expression  and  necessary  references. 

Appendix  C provides  a physical  description  of  the  underwater  vehicle 
used  as  the  basis  for  calculations  illustrated  graphically  in  this  report, 
and  Appendix  D contains  time  history  plots  of  the  four  maneuvers  from 
which  the  trajectory  simulation  and  percent  contribution  plots  were  ob- 
tained. Reference  frame  conventions  are  as  shown  in  Figure  1. 

A word  of  caution  is  appropriate  here.  When  comparing  computed  ac- 
celeration coefficient  values  with  experimentally  determined  values, 
attention  must  be  given  to  the  location  of  the  reference  frame  origin. 

If  the  origin  during  the  experimental  tests  is  not  the  vehicle  eg,  then 
erroneous  comparisons  will  result  for  those  coefficients  where  distances 
from  the  eg  are  involved. 


LONGITUDINAL  COEFFICIENTS 


1.  X’ . is  an  acceleration  coefficient  resulting  from  the  resisting 
force  of  the  fluid  on  a vehicle  accelerating  in  the  x-direction.  This 
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force  is  required  to  change  the  flow  field  about  the  vehicle,  whereas 
the  static  coefficients  represent  the  forces  necessary  to  maintain  a 
flow  field  configuration  at  a constant  vehicle  velocity.  More  simply, 
X'^  may  be  thought  of  as  a mass  of  fluid  that  must  be  accelerated  along 
with  the  vehicle,  thus  giving  the  vehicle  an  "added"  or  "apparent"  mass 
which  is  manifest  during  vehicle  acceleration. 

2.  Insight  to  the  relative  significance  of  X’.  can  be  gained  by 
noting  that  X' ^ appears  in  the  linear  equations  of  motion  in  the  mass 
term 


(m*  - X’.)  . 
u 

Thus,  in  cases  where  the  value  of  X*.  represents  a significant  contribu- 
tion to  this  term,  accuracy  in  predicting  X'.  is  important.  The  X'. 
contribution  to  vehicular  moti  ' during  a control  surface  step  input 
is  generally  minor,  and  dissipai.es  rapidly  as  the  vehicle  attains  a 
steady  state.  In  control  surface  reversal  maneuvers,  as  Illustrated 
in  trajectory  simulation  and  x-force  percent  contribution  plots  in 
Figures  2 and  3,  X’ . becomes  slightly  more  significant  and  longer-acting, 
but  still  remains  a relatively  minor  influence.  For  nearly  neutrally 
buoyant  faired  vehicles  having  an  £,/ d ratio  higher  than  8.0,  X’.  typi- 
cally contributes  less  than  3 percent  to  the  mass  term.  X*.  has  in- 
creasing relative  significance  for  lower  £/d  ratios.  Root  iocus  studies 

for  typical  submarine  and  torpedo-like  bodies  has  corraborated  the  fore- 
going, indicating  very  minor  sensitivity  to  large  percentage  errors  in 
predicting  X'^. 

When  performing  a two-degree-of-freedom  linear  analysis  in  the  lon- 
gitudinal plane  (M  and  Z equations),  three  roots  of  the  characteristic 
equation  are  obtained.  The  investigator  may  find  a fourth  root  arising 

from  the  x-degree  of  freedom  equation  by  assuming  that  the  only  signifi- 

cant terms  are  X'u  and  X\  (this  is  often  acceptable  for  an  axisymmetric 
vehicle) . The  fourth  root  is  found  by  the  relation 


X' 


Tables  1A  through  ID  present  the  results  of  a sensitivity  analysis  of 
nondimensional  roots  of  the  characteristic  equation  (transfer  function 
denominator)  and  u,  w,  and  9 (numerators)  for  a generalized  underwater 
vehicle,  varying  X\  from  -100  percent  to  +100  percent. 

3.  Mathematically,  given  the  linearizing  truncation  of  the  Taylor 
series  expansion  to  the  first  term,  X\u  is  the  force  additional  to  the 
inertial  force  occurring  with  a vehicular  x-acceleration.  In  the  com- 
putation of  X\,  the  body  is  approximated  by  a prolate  ellipsoid  having 
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like  length  and  maximum  diameter,  and  neglecting  appendages  to  the  body 
due  to  their  comparatively  small  frontal  area.  Using  the  results  of 
Lamb^1^,  derived  assuming  potential  flow,  we  have: 


4. 


kimdf 

W ’ 


(11) 


where 


with 

a0  " <^n  [l4i]  ” e)  ’ (13) 

e being  the  eccentricity  of  the  rotated  ellipse  expressed  in  common 
terms  as 


kl  " 


2 - a. 


(12) 


e 


(14) 


The  inertia  coefficient  for  axial  flow,  kj_,  ranges  from  0.5  to  0 
for  fineness  ratios  of  1.0  to  09 . k.^  is  graphically  presented  in  Figure  4. 

Z'  . 

u 

1.  The  acceleration  coefficient  Z'.  represents  a z-force  resulting 
from  an  acceleration  in  the  x-direction. 

2.  Root  locus  sensitivity  studies  on  typical  submarine  and  torpedo- 
like vehicles  show  extreme  insensitivity  to  large  magnitude  excursions 
and  sign  changes  of  Z'.  . Thus  the  coefficient  is  generally  assumed 
negligible  and  set  to  zero. 

3.  Vehicles  exhibiting  symmetry  about  the  xy-plane  produce  balanced 
z-forces  during  acceleration  along  x,  indicating  Z'.  - 0.  In  potential 
theory,  symmetry  about  the  xy-plane  or  yz-plane  causes  Z*.  ■ 0 b.s  shown 
in  Appendix  A. 

4.  In  the  vast  majority  of  cases,  the  combination  of  near  symmetry 
about  the  xy-plane  and  minor  influence  on  roots  of  the  characteristic 
equation  leads  to  setting  Z’^  *»  0. 


(l)ibid. 


(Text  Continued  on  Page  16) 
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FIGURE  4.  LAMB'S  INERTIA  COE-FICIENTS  k,,  AND  k'  PLOTTED  AS  A FUNCTION  OF  VEHICLE  t/d 
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1.  The  term  M' . represents  the  pitching  moment  resulting  from  an 
acceleration  parallel  to  the  x-axis. 


2.  Root  locus  sensitivity  studies  on  typical  submarine  and  torpedo- 
like bodies  show  extreme  insensitivity  to  large  magnitude  changes  in  M\  . 

3.  As  was  the  case  with  Z'.,  vehicles  having  symmetry  or  near  symme- 
try about  the  xy-plane  will  create  a negligible  M-moment  due  to  x- 
acceleration.  Potential  theory  predicts  ■ 0 for  yz-plane  symmetry 
(see  Appendix  A).  Extreme  vertical  separation  of  the  cb  and  eg,  coupled 
with  a large  X'.  force  (arising  from  low  H/d  ratio),  could  conceivably 
cause  a significantly  large  M' . . For  free-flooded  vehicles,  launcher  de- 
vices, and  torpedo-like  bodies,  this  is  a highly  unlikely  occurrence. 


4.  The  combination  of  near  symmetry  about  the  xy-plane  and  minor 

influence  on  the  roots  of  the  characteristic  equation  usually  leads  to 

setting  Mr . - 0.  For  a vehicle  with  an  extremely  large  vertical  eg  - cb 

separation  and/or  an  X' . value  that  is  significant  when  compared  to  m' , 

M’ . is  calculated  as 
u 


M’ 


« 

u 


< 


where  ■ z-distance  from  the  eg  to  the  cb,  positive  down. 


X' 


w 


1.  The  acceleration  derivative  X'^  represents  the  x-force  resulting 
from  an  acceleration  in  the  z-direction. 

2.  As  with  Z* . and  M’ . , the  roots  of  the  characteristic  equation 
show  negligible  change  with1^' . set  equal  to  zero,  then  varied’  positively 
and  negatively  over  a range  much  greater  than  would  be  likely  in  reality. 

3.  Near  balanced  z-force  during  z-acceleration  due  to  near  fore- 
aft  symmetry  yields  resultants  of  a large  z-force  component  and  some 
pitching  moment,  but  very  little  x-force.  It  is  significant  to  note 
that  Z' . «*  X’ . according  to  potential  theory  (Equation  (9)).  If  the  eg 
and  cb  are  significantly  separated  along  both  the  x-  and  z-axes  then, 
since  those  forces  arising  from  potential  flow  theory  are  in  a reference 
frame  with  a cb  origin,  the  pitching  moment  term  M'^  will  yield  an  X'^  term 
when  transferred  to  the  eg.  As  with  Z^,  though,  even  variation  through 
extreme  ranges  showed  that  setting  X' . ■ 0 was  justified  for  application 

to  the  vast  majority  of  underwater  vehicles. 


4. 


0. 
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1.  The  resisting  hydrodynamic  force  resulting  from  a *- 
acceleration  is  represented  in  the  equation  of  motion  by  Z'^  . The 
term  (m'  - Z'  .)  is  the  nondimensional  "apparent  mass"  of  the  vehicle 
when  one  is  considering  heavn  acceleration. 


2.  Z'^  is  one  of  the  most  important  acceleration  derivatives. 
Trajectory  simulation  and  z-force  percent  contribution  plots  in  Fig- 
ures 5 through  10  illustrate  its  significance  to  vehicular  motion 
during  various  maneuvers.  The  root  locus  plots  in  Figures  11,  12,  and 
13  illustrate  the  wide  variation  in  the  locations  of  the  poles  due  to 
variation  in  Z'..  Tables  2A  through  2D  present  the  results  of  a sensi- 
tivity analysis*}.*  nondimensional  roots  of  the  characteristic  equation 
(transfer  function  denominator)  and  u,  w,  and  0 (numerators)  for  a gen- 
eralized underwater  vehicle,  varying  Z'.  from  -100  percent  to  +100  per- 
cent. As  in  X* . and  M' . , insight  into  \:he  influence  of  Z'^  can  be  gained 
by  noting  the  significance  of  Z'.  in  the  hydrodynamic  mass  term  (m*  - 
Z'.).  Vehicles  with  i/d  ratios  greater  than  8.5  will  usually  have 
nearly  a 50  percent  Z'.  contribution  to  the  apparent  mass.  That  is,  the 
added  mass  will  be  roughly  equal  to  the  mass.  Horizontal  fins  will  con- 
tribute to  Z'.  and  frequently  cause  it  to  be  larger  than  m* . 

* ’’’  l 

3.  As  with  $ . ^.he  body  is  approximated  by  a prolate  ellipsoid 
having  like  lengtn  an#  maximum  diameter.  Lamb's^15  inertial  coefficient 
for  such  a body  in  ctyisswise  flow  is  used  as  follows: 


Z' . 
wv 


where 


k2mdf 

*SPi3 


k2V 

Iji3 


2 - 6q  ’ 


(15) 


(16) 


with 


- -tt  - 


1 - e 
2e3 


in 


1 + e 
1 - e 


(17) 


e being  the  eccentricity  of  the  rotated  ellipse  as  noted  in  Equation  (14) . 
k2  is  graphically  presented  in  Figure  14.  Fortunately,  this  very  impor- 
tant body  coefficient  has  proved  to  be  accurate  to  ±5  percent  in  compari- 
son with  experimental  data  on  the  majority  of  underwater  vehicles 


( 1 ) ibid. 


(Text  Continued  on  Page  23) 
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FIGURE  7.  TRAJECTORY  SIMULATION  PLOT  OF  Z'.  CONTRIBUTION  TO  Z-FORCE 
DURING  A «$  - ±30*  DIVE/CLIMB  MANEUVER 
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FIGURE  8.  PERCENT  CONTRIBUTION  OF  Z'^  TO  Z-FORCE  DURING  A 
«.  ■ ±30*  D I VC/CL  I Mi  MANEUVER 
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compared.  In  application  to  a wide  variety  of  underwater  vehicles,  it 
has  been  our  experience  that  if  a body  is  nearly  axisynmetric  then 
Lamb's  inertial  coefficient  is  a reliable  formulation  for  Z'.  . This 
is  especially  true  for  large  i/d  ratios.  b 

This  is  the  first  hydrodynamic  coefficient  with  significant  con- 
tribution from  appendages  to  the  body.  The  discussion  and  formulations 
that  follow  can  be  applied  to  any  surface  projecting  from  the  body  that 
can  be  approximated  as  a flat  plate  in  accelerated  fluid  flow  normal  to 
its  surface.  The  equations  given  will  be  for  horizontal  tail  fins,  but 
they  are  applicable  to  any  horizontal  surfaces  such  as  bowp lanes,  sail- 
planes, fairwater  planes,  depressors,  wings,  faired  struts,  sails,  etc. 
Tail  contributions  to  other  longitudinal  hydrodynamic  coefficients  are 
formulated  using  the  expressions  given  below.  Theoretical  values  of  the 
additional  mass  of  a number  of  types  of  bodies  (bodies  of  infinite  length, 
bodies  of  revolution,  ellipsoids  or  elliptic  plates  of  finite  dimensions) 
have  been  derived  by  Lamb^1^  and  Munk^2^3^.  Values  for  bodies  of  finite 
dimensions  not  covered  by  the  theory;  for  example,  rectangular  plates  and 
wings,  have  been  the  object  of  many  experimental  research  programs  on  the 
phenomenon. 

Extensive  test  programs  on  additional-mass  and  additional-moment 
effects  were  conducted  by  the  National  Advisory  Committee  for  Aero- 
nautics (NACA)  from  1933  to  1940.  The  coefficients  given  here  are  based 
on  the  results  of  these  tests  as  reported  by  Gracey^4 * * *^ , and  Malvestuto 
and  Gale(B) . A very  similar  formulation  is  given  by  Pastor  and 
Abkowitz^8) . 

For  a thin  flat  plate  of  infinite  span  moving  in  a perfect  fluid  at 
constant  velocity,  V,  along  the  normal  to  its  surface,  the  momentum  im- 
parted to  the  fluid  per  unit  span  is  given  by  hydrodynamic  theory  as 

2„ 

irpc  V 

4 ‘ 

For  plates  of  finite  span,  this  expression  must  be  corrected  by  the 
introduction  of  coefficients  whose  values  depend  on  the  dimensions  of 


(1)ibid. 

(a)ibid. 

(a)ibid. 

<4)ibid. 

<B>ibid. 

^Naval  Underwater  Ordnance  Station  T.N.  No.  120,  Hydrodynamic  Stability 

and  Control  Derivatives,  by  D.  L.  Pastor  and  M.  A.  Abkowitz,  1957. 
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the  plate.  The  additional  mass,  m , for  translation  of  a plate  of  span 
b,  V is  thus  determined  from  the  equation  of  linear  momentum 

„ _ kTrpc2bV 

tn  v — z ) 

a 4 f 


so  that 


m 

a 


2 

kirpc  b 
4 


(18) 


where  k is  the  coefficient  of  additional  mass  shown  in  Figure  14.  This 
formulation  is  for  one  fin.  So,  with  the  proper  nondimensionalizing 
terms  added,  the  expression  for  two  identical  horizontal  fin  surfaces 
in  a tail  configuration  is  written 


2 2 

1 kirpc  b kirc  b 


(19) 


Note:  Pastor  and  Abkowitz  substituted  a curve  fitted  expression  for  k: 


(1  + l/AR)^ 


An  effective  chord,  c,  is  usually  used  since  most  surfaces  are  not  simple 
rectangular  shapes.  If  a sternplane  or  any  movable  horizontal  tail  sur- 
face is  present  and  is  an  extension  of  the  fixed  horizontal  tail  surface, 
then  it  should  be  included  in  the  area  for  calculating  c and  s.  If  it  is 
not  a part  of  the  fin,  then  its  contribution  to  Z'.  can  be  calculated 
separately,  as  can  the  bowplane  contribution,  by  using  Equation  19  to 
yield  Z' . and  Z' . . For  completeness  and  clarification.  Figure  15  and 

s bp 

the  following  relationships  are  given. 

2 

AR  = aspect  ratio  of  surface  = b /s 


b = root-to-tip  span  of  surface 


s = planform  area  of  surface 


c = root  chord  of  surface 
r 

cfc  = tip  chord  of  surface 
c = mean  chord  of  surface  = s/b  = b/AR 


= planform  taper  ratio  of  surface  = c /cr 


(Text  Continued  on  Page  26) 
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Tables  3A  through  3D  present  the  results  of  a sensitivity  analysis  of 
nondimensional  roots  of  the  characteristic  equation  (transfer  function 
denominator)  and  u,  w,  and  6 (numerators)  for  a generalized  underwater 
vehicle,  varying  M' . from  -100  percent  to  +100  percent. 


3.  The  eg  and  cb  of  a submerged  vehicle  are  generally  in  vertical 
alignment  to  i chieve  level  trim  at  zero  speed.  Since  the  body  addi- 
tional mass  effects  are  lumped  as  acting  through  the  body  cb  (caution 
must  be  exercised  when  the  cb  is  not  coincident  with  the  centroid  of 
the  volume  displaced  by  the  outer  hull  surface — as  with  "wet"  vehicles) 
and  the  vehicle  equations  of  motion  are  typically  written  about  a eg 
origin,  the  resisting  force  to  a z-acceleration  is  a pure  Z'.  w term. 

Wb 


That  is,  there  is  no  x-moment  arm  to  cause  the  generation  of  a moment 
r*  ” The  added  mass  and  moment  effects  of  the  horizontal  surfaces 


M'  . w. 
w 

generally  are  tha  primary  contributors  to  M* 


moment 

arm  is  calculated  to  produce 

4. 

M\ 

1 

- z' . - 

wb 

Wb 

M’ . 

- Z*.  - 

w 

t 

ft 

w 


Just  as  any  force  on 
30  with  M' . we  have 


(X  . - X ) 
neb  neg 


(Xnht  - W 


(21) 


(22) 


where 


neb 

= x-distance 

from 

nose 

to 

cb 

neg 

■ x-distance 

from 

nose 

to 

eg 

nht 

= x-distance 

from 

nose 

to 

centroid 

of  the  horizontal  tail 

fins  . 

Note:  The  distances  are  all  positive  numbers. 


Substituting  the  x-distance  to  bowplanes  or  shroud  into  Equation  (22) 
yields  M'^  and  M'^  . The  complete  expression  is  the  sum  of  the  com- 

ponent contributions f e.g., 


M' 


= M' 


w. 


M\ 

w, 


ht 


M\ 

w, 


bp 


If. 

w 
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1.  Forces  along  the  x-direction  caused  by  angular  acceleration  in 
pitch  are  represented  by  X\  . 

2.  & 3.  X'^  is  usually  an  insignificant  term  as  is  its  near-equal, 

M\  (Equation  (9)). 


4.  X' . = 0 . 

q 


1.  z-forces  arising  from  pitch  angle  accelerations  are  termed  Z'. 

q 

2.  Z'.  , as  its  counterpart  M'^  , is  a minor  derivative,  yet  it 

can  have  soSe  influence  over  the  characteristic  equation  of  a vehicle. 
It  should  therefore  be  calculated.  Tables  4A  through  4D  present  the 
results  of  a sensitivity  analysis  of  nondimens ional  roots  of  the  char- 
acteristic equation  (transfer  function  denominator)  and  u,  w,  and  0 
(numerators)  for  a generalized  underwater  vehicle,  varying  Z'.  from 
-100  percent  fo  +100  percent.  ^ 

3.  & 4.  Z\  = M'  . . 

q W 


1.  When  a submerged  body  is  accelerated  in  pitch,  it  must  change 
the  flow  pattern  in  the  surrounding  fluid.  The  added  moment  needed  to 
effect  this  change  in  flow  pattern  is  represented  in  the  equations  of 
motion  as  an  apparent  moment  of  inertia,  (I'  - M'.)  . 


2. 


M' . is  a major  derivative  in  influence  on  the  characteristic 


equation.  As  with  Z* . and  X'.,  M'  , enters  the  equations  of  motion  as 
an  addition  to  the  vehicle  inertia^term,  in  this  case  the  moment  of 
inertia  about  the  y-axis.  Trajectory  simulation  and  pitch  moment 
percent  contribution  plots  in  Figures  16  through  21  illustrate  its  sig- 
nificance to  vehicular  motion  during  various  maneuvers.  The  root  locus 
plots  in  Figures  22,  23,  and  24  illustrate  the  wide  variation  in  the 
locations  of  the  poles  due  to  variations  in  M' . . laoles  5A  through  5D 


present  the  results  of  a sensitivity  analysis  of  nondimens ional  roots 
of  the  characteristic  equation  (transfer  function  denominator)  and  u, 
w,  and  8 (numerators)  for  a generalized  underwater  vehicle,  varying  M 
from  -100  percent  to  +100  percent. 


q 


(Text  Continued  on  Page  36) 
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3.  & 4.  The  body  contribution  to  M\  is  calculated  assuming  the 

body  is  a prolate  ellipsoid.  Thus,  we  use  Lamb's  coefficient  for 
added  moment  of  inertia  and  the  parallel  axis  theorem 


-k'  1 
b y 


M' 


df 


ho 


+ z' 


(X  . - X )' 
neb  neg 

F 


where 


k'  = 
K b 


e (3q  - aQ) 


(2  - e2)  {2e2  - (2  - e2)(BQ  - aQ)} 


aQ  and  6q  are  defined  in  Equations  (13)  and  (17).  These  equations  are 
repeated  here  for  convenience: 


2(1  - e2)  „ „ fl  + el 

a0  - — 3 — (is£n  Lr^J ' e)  ’ 


g _ 1_  _ 1.  - e £n  + e 

30  ~ 2 0 3 A 1 - e ’ 

e 2e 


(13) 

(17) 


[ = moment  of  inertia  about  the  y-axis  of  the 

ydf  mass  of  the  fluid  displaced  by  the  body. 


The  tail  contribution  to  M' . is  in  the  same  form  as  the  body  terms.  That 

q 

is,  there  is  an  added  moment  of  inertia  term  and  a term  representing  an 
added  mass  acting  on  a moment  arm  from  the  eg.  The  added  moment  of  in- 
ertia of  a flat  plate  resulting  from  rotation  about  an  axis  in  the  plane 
of  the  plate  and  parallel  to  the  span  is  given  by 


I = 
a 


i . v2“3 

k'  irpb  c 

_£ 

48 


where  k'  is  the  coefficient  of  added  moment  of  inertia  for  a flat  plate 
plottedpon  Figure  25.  When  c/b  is  sufficiently  small  (AR  sufficiently 
large),  this  term  becomes  negligible  compared  with  the  second  term  in 
M' . given  as 

"t 

z,.  (iInht  - W 


W, 


ht 


V 
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The  horizontal  tail  surface  contribution  is  then 


M'  . 


. (k'  irpb2c3) 

1 X 


*ht 


p45  48 


(X  , - X )‘ 

+ z\  ■ — .n(=s. 


w, 


ht 


Similar  expressions  may  be  written  for  bowplanes  and  a shroud  so  that 
:omplete  M\  formula  would  be 


M'  . = M'.  + M’ . + M'.  + M* 


• • 


Tit 


bp 


*sh 


LATERAL  COEFFICIENTS 


Y' 


v 


1.  Forces  along  the  y-axis  proportional  to  accelerations  along 
the  same  axis  are  characterized  by  the  added  mass  term  Y'.  . This  term 
is  the  lateral  counterpart  of  Z'^  in  the  longitudinal  equations. 

2.  Just  as  Z'.  was  seen  to  be  a major  stability  derivative  in  the 
longitudinal  plane, WY*.  is  a very  important  coefficient  in  the  lateral 
characteristic  equation.  Trajectory  simulation  and  y-force  percent  con- 
tribution plots  in  Figures  26  through  29  illustrate  its  significance  to 
vehicular  motion  during  various  maneuvers.  The  root  locus  plots  in  Fig- 
ures 30,  31,  and  32  illustrate  the  wide  variation  in  the  locations  of 
the  poles  due  to  variation  in  Y\  . Tables  6A  through  6D  present  the 
results  of  a sensitivity  analysis  of  nondimensional  roots  of  the  char- 
acteristic equation  (transfer  function  denominator)  and  v,  ip,  and  <|> 
(numerators)  for  a generalized  underwater  vehicle,  varying  Y'^  from 
-100  percent  to  +100  percent. 

3.  & 4.  For  axisymmetric  bodies 


Most  underwater  vehicles  have  nearly  axisymmetric  hull  forms.  Of  those 
whose  shapes  are  not  a body  of  revolution,  quite  often  good  approxima- 
tions can  be  made  by  treating  the  body  as  everywhere  circular  in  cross- 
section  with  a cross-sectional  area  equal  to  the  actual  body.  This 
approximation  is  often  used  in  the  analysis  of  underwater  vehicles. 


Appendages  that  contribute  to  Y\  are  those  with  significant  pro- 
jections onto  the  xz-plane.  Usually  these  consist  of  vertical  tail  fins 


(Text  Continued  on  Page  43) 
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(or  stabilizers),  rudders,  canards,  and  shrouds.  The  contributions  of 
these  and  any  other  surfaces  that  can  be  treated  as  flat  plates  in 
accelerated  fluid  flow  normal  to  its  surface  can  be  calculated  using 
Equation  (18).  It  is  repeated  here  for  convenience. 


m 

a 


knpc  b 


(18) 


Since  it  is  common  to  have  unequal  upper  and  lower  vertical  fins 
in  the  tail,  the  two  will  not  be  combined  as  were  the  two  horizontal 
tail  fins.  The  expressions  for  the  upper  and  lower  vertical  tail  fins 
are  written  as 


Y' . 

vuvt(ivt) 


kirc^b 


2T 


If  the  rudder  is  an  extension  of  the  vertical_tail  fins,  then  it 
should  be  included  in  the  area  for  calculation  of  c.  If  the  rudder  is 
separate,  them  it  can  be  treated  separately  using  the  same  expression. 

Since  the  bowplanes  have  no  significant  vertical  surface  area, 
they  do  not  contribute  to  Y' . . The  shroud  has  the  same  projection  in 
the  vertical  plane  as  in  thevhorizontal  plane,  30 


The  complete  expression  of  Y'^  for  a typical  underwater  vehicle  is 
written 


+ Y* . + Y' . 

V V 

uvt  ivt 


+ 


1.  y-forces  on  a submerged  vehicle  that  are  due  to  angular  accel- 
eration in  roll  are  represented  in  the  equations  of  motion  as  Y'.  . 

2.  This  term  can  be  significant  on  vehicles  with  large  Z^/l, 

Z /£,  and  Z.  Ji.  ratios,  where 
uct  £ct 

Z , - z-distance  from  body  eg  to  cb 

cb 

Z _ - z-distance  from  the  x-axis  to  the  centroid 
uct 

of  the  upper  vertical  tail,  and 

Z.  ■ z-distance  from  the  x-axis  to  the  centroid 
c of  the  lower  vertical  tail. 
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Tables  7A  through  7D  present  the  results  of  a sensitivity  analysis  of 
nondimensional  roots  of  the  characteristic  equation  (transfer  function 
denominator)  and  v,  ij>,  and  <j>  (numerators)  for  a generalized  underwater 
vehicle,  varying  Y\  from  -100  percent  to  +100  percent. 

3.  Due  to  intermodal  coupling,  Y'.  - K' . as  shown  in  Equation  (8). 
Since  the  conceptual  description  of  the^physical  phenomenon  causing  K' . 
moments  is  more  easily  understood  than  that  for  Y' . , only  the  equality 
relationship  is  given  here;  that  is  P 

Y' . ■ K' . (See  below), 
p v 


When  an  angular  acceleration  in  yaw  results  in  a force  in  the 
y-direction,  the  y-force  is  represented  in  the  equations  of  motion  by 
the  term  Y'^  . The  longitudinal  counterpart  of  Y\  is  Z\  . 

Tables  8A  through  8D  present  the  results  of  a sensitivity 
analysis  of  nondimensional  roots  of  the  characteristic  equation  (trans- 
fer function  denominator)  and  v,  <Ji,  and  <j>  (numerators)  for  a generalized 
underwater  vehicle,  varying  Y* . from  -100  percent  to  +100  percent. 


3.  Intermodal  coupling  yields  the  relationships  in  Equation  (8). 
That  is. 


4. 


1.  Rolling  moment  due  to  acceleration  along  the  y-axis  is  termed 


2.  K'  . couples  the  roll  and  yaw  motions  of  a vehicle.  As  a 
vehicle  enters  a turn,  the  changes  in  the  drift  angle,  8,  will  appear  as 
a change  in  the  sway  velocity,  v;  i.e.,  an  acceleration  v.  The  rolling 
moment  produced  by  v in  this  situation  can  be  significant  on  a vehicle 
with  unsymmetrical  vertical  fin  arrangements,  like  a submarine's  sail. 
Trajectory  simulation  and  roll  moment  percent  contribution  plots  in 
Figures  33  and  34  illustrate  its  significance  to  vehicular  motion  during 
a turn  reversal  maneuver.  The  root  locus  plot  in  Figure  35  illustrates 
the  variation  in  the  locations  of  the  poles  due  to  variation  in  K'^  . 
Tables  9A  through  9D  present  the  results  of  a sensitivity  analysis  of 
nondimensional  roots  of  the  characteristic  equation  (transfer  function 
denominator)  and  v,  ijj,  and  <j>  (numerators)  for  a generalized  underwater 
vehicle,  varying  K\  from  -100  percent  to  +100  percent. 


(Text  Continued  on  Page  50) 
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FIGURE  35.  ROOT  LOCUS  PLOT  OF  4,/S  FOR  A GENERALIZED  UNDERWATER  VEHICLE,  VARYING  K’ 
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3.  & 4.  This  coefficient  is  evaluated  in  a similar  fashion  to 
that  of  M*  . . It  is  the  representation  o£  added  mass  effects  acting 
on  moment  arms  from  the  eg. 


1 

As  noted  before,  it  is  not  uncommon  for  the  upper*  and  lower  verti- 
cal tail  fins  to  be  of  different  size  and  shape.  Also',  since  the  eg  is 
generally  below  the  body  axis  of  symmetry,  there  are  unequal  moment 
arms  to  the  centroids  of  the  upper  and  lower  fins  even  if  the  fins  are 
symmetric  about  the  center  line.  The  vertical  tail  surfaces’  contri- 
bution to  K’ . is  written 
v 


uct 


+ Y’. 

v 

£vt 


z 


ct 


} 


= K’ . + K’  . 

V V 

uvt  £vt 


Z ^ distance  from  the  x-axis  to  the  centroid 
uct 

of  the  upper  vertical  tail 

= distance  from  the  x-axis  to  the  centroid 
c of  the  lower  vertical  tall. 


Note:  The  vertical  tail  fins  are  assumed  to  lie  in  the  xz-plane. 


The  complete  expression,  then,  is  written 


K’ 


v 


Any  other  vertical  surfaces  that  can  be  approximated  as  flat  plates 
lying  in  the  xz-plane  can  be  treated  similarly  to  the  upper  and  lower 
vertical  tail  surfaces. 


K\ 


1.  Moments  additional  to  the  Ix  moment  of  inertia  resisting  angu- 
lar acceleration  in  roll  are  included  in  the  term  K’ . . This  added 
moment  of  inertia  is  one  that  has  no  Lamb  coefficient1^  for  the  body 
contribution.  This  is  because  a homogeneous  body  of  revolution  can 
rotate  about  its  axis  of  symmetry  and  not  experience  any  forces  from  an 
inviscid  fluid. 


(Dibid. 
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2.  Tables  10A  through  10D  present  the  results  of  a sensitivity 
analysis  of  nondimens ional  roots  of  the  characteristic  equation  (trans- 
fer function  denominator)  and  v,  and  <j>  (numerators)  for  a generalized 
underwater  vehicle,  varying  K'^  from  -100  percent  to  +100  percent. 

3.  Since  p,  angular  acceleration  in  roll,  is  defined  about  the 
x-axis  passing  through  the  eg,  there  will  be  some  added  moment  effect 

of  the  body  due  to  an  added  mass  effect  acting  at  the  cb.  It  is  written 
as 


Any  surface  projecting  radially  from  the  vehicle  (such  as  tail  fins, 
bowplanes,  etc.)  will  contribute  to  the  added  moment  effect  represented 
by  K' . . Each  surface  contributes  in  two  ways.  One  is  as  an  added  mass 
acting  on  a moment  arm  from  the  x-axis.  The  general  form  for  this  is 
the  same  as  the  body  contribution  above.  Considering  any  "flat  plate" 
surface,  such  as  the  upper  vertical  tail  fin,  this  term  would  be  written 


K'. 


= Y'. 


uct 


uvt 


uvt 


The  other  form  of  contribution  to  K' . is  by  the  individual  added  moment 
of  inertia  of  the  flat  plate  rotating  about  an  axis  parallel  to  the 
x-axis  but  passing  through  the  centroid  of  the  flat  plate.  Again,  we 
turn  to  Malvestuto  and  Gale^6^  for  experimental  data  on  flat  plates  with 
varying  dihedral  angle  and  taper  ratio.  For  rotation  about  an  axis  in 
the  plane  of  the  plate,  parallel  to  the  chord,  and  passing  through  the 
centroid,  the  equation  is 


I 

x 

a 


TIP 


2,  3 


D,D„k'cb 


48  A*T 


where  k'  is  the  coefficient  of  added  moment  of  inertia  for  a flat  plate 
plotted*5 in  Figure  25,  is  the  correction  factor  for  taper  ratio  plotted 
in  Figure  36,  and  Dp  is  the  correction  factor  for  dihedral  angle  plotted 
in  Figure  37.  Thus,  for  any  "flat  plate"  surface  such  as  the  upper  verti- 
cal tail  fin 


K' 


uvt 


*5pJl5  Ix 


+ Y’ 


uct 

nr 


(6)ibid. 
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FIGURE  37.  VARIATION  OF  THE  ADDITIONAL  MOMENTS  OF  INERTIA 
OF  A SINGLE  PLATE  WITH  DIHEDRAL  ANGLE 


) 

E 

I 


The  first  term  is  often  negligible  compared  to  the  second.  The  complete 
expression  for  a common  vehicle  would  be 


1.  Rolling  moments  caused  by  angular  accelerations  in  yaw  are  rep- 
resented by  K'.  . 

r 

2.  Tables  11A  through  11D  present  the  results  of  a sensitivity  analy- 
sis of  nondimens ional  roots  of  the  characteristic  equation  (transfer  func- 
tion denominator)  and  v,  iji,  and  4>  (numerators)  for  a generalized  underwater 
vehicle,  varying  K' ^ from  -100  percent  to  +100  percent. 

3.  K*  . arises  in  the  same  manner  as  the  K\v  term  except  that  th ; 
linear  acceleration  normal  to  the  vertical  surfaces  is  the  result  of  an 
angular  acceleration  in  yaw  acting  about  the  z-axis. 
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K',  - K*.  ^nuvt  ~ Xncg^  K,  ^Xnftvt  ~ Xncg^ 

1 ♦»«  ‘ 

where 

X^uvt  “ distance  from  nose  to  upper  vertical 
tail  surface  centroid,  positive  aft 

X * ^^-stance  fro®  nose  to  lower  vertical 
tail  surface  centroid,  positive  aft. 


1.  Linear  accelerations  in  the  y-direction  causing  yawing  moments 
give  rise  to  the  stability  derivative  N'.  to  describe  these  resulting 
moments . 


2 .  Tables  12a  through  12D  present  the  results  of  a sensitivity 
analysis  of  nondimensional  roots  of  the  characteristic  equation  (trans- 
fer function  denominator)  and  v,  , and  <j>  (numerators)  for  a generalized 
underwater  vehicle,  varying  N*  . from  -100  percent  to  +100  percent. 


3.  Primarily  due  to  the  vertical  tail  surfaces,  a vehicle  changing 
drift  angle  will  experience  a restoring  force  proportional  to  v (result- 
ing from  &).  N\  is  the  term  representing  this  force. 


4. 


(X  . - X ) 
neb  neg 

l 


N'  . 

v. 


-Y\ 

v 


(X  - X ) 
nuvt  neg7  _ , 


uvt 


ivt 


1. 

resented 


Yawing  moments  caused  by  angular  accelerations  in  roll  are  rep- 
in  the  vehicular  equation  of  motion  by  N \ . 


2.  Tables  13A  through  13D  present  the  results  of  a sensitivity 
analysis  of  nondimensional  roots  of  the  characteristic  equation  (trans- 
fer function  denominator)  and  v,  4>,  and  $ (numerators)  for  a generalized 
underwater  vehicle,  varying  from  -100  percent  to  +100  percent. 

3.  & 4.  Equation  (9)  show3 
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1.  Yawing  moments  additional  to  I are  denoted  N*.  . This  is  the 
added  moment  of  inertia  term  whose  longitudinal  counterpart  is  . 

2.  Trajectory  simulation  and  yaw  moment  percent  contribution  plots 
in  Figures  38  and  39  illustrate  its  significance  to  vehicular  motion 
during  various  maneuvers.  The  root  locus  plots  in  Figures  AO,  41  and 

42  illustrate  the  wide  variation  in  the  locations  of  the  poles  due  to 
variation  in  N*.  . Tables  14A  through  14D  present  the  results  of  a sen- 
sitivity analysis  of  nondimensional  roots  of  the  characteristic  equation 
(transfer  function  denominator)  and  v,  ip , and  $ (numerators)  for  a gen- 
eralized underwater  vehicle  varying  N* . from  -100  percent  to  +100 
percent.  r 


3.  & 4.  Like  M*.  , there  are  two  parts  to  the  N*.  body 
contribution  ” r 


KbIz,f  (X  - X )' 

N\  - — i P-  + Y'  - ncb... "c8 

rb  Y vb  l2 


for  an  axisymmetric  vehicle  body  1 - I . Again,  similar  to  M'.  , 


“df 


df 


the  upper  vertical  tail  fin  contribution  is 


uvt 


(K^TTPb2C3) 

"W  48 


+ Y'. 


^Xnuvt-  Xnc/ 

V . I2 

uvt 


The  lower  vertical  tail  fin  is  identical  in  form.  The  complete  expres- 
sion, then,  is  written  N?.  = N*.  + II' . + N'. 

1 rb  ruvt  r£vt 


SUMMARY 


A compilation  of  methods  for  analytically  predicting  the  accelera- 
tion hydrodynamic  coefficients  from  the  geometric  and  mass  distribution 
characteristics  of  underwater  vehicles  has  been  presented.  A theoreti- 
cal development  from  potential  flow  theory  is  given  to  show  the  equiva- 
lences among  acceleration  and  coupled  terms.  The  detailed  discussion  of 
each  coefficient  presented  includes  descriptions  intended  to  provide  the 
reader  with  an  understanding  of  the  physical  relationships  involved. 
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Results  of  sensitivity  analyses  are  presented  to  show  the  relative  con- 
tributions of  each  coefficient  to  the  equations  of  motion.  A discussion 
of  the  important  assumptions  involved  in  mathematically  computing  each 
coefficient  is  given,  along  with  the  analytical  expression  used  to  com- 
pute each  coefficient.  An  example  case  (Appendix  B)  is  given  to  show  a 
comparison  between  computed  coefficients  and  experimental  data. 

Table  13  presents  a summary  of  the  sensitivity  analyses  conducted 
for  each  acceleration  coefficient.  These  results  are  for  ±100  percent 
coefficient  variation.  It  can  be  seen  in  Table  15  that  five  coefficients 
make  major  contributions  to  the  vehicle  dynamics.  These  are  Z'^,  M'., 

Y' . , K' . , and  1 . . Variations  of  ±100  percent  change  both  the  numerator 
anti  denominator  roots  by  better  than  50  percent,  except  for  K'^  . Changes 

in  K* . significantly  effect  only  the  roll  numerator,  It  should  be 

v r 

noted  that  the  v terms,  such  as  K' . and  Y'.,  do  not  affect  the  v-numerator 
because  no  v-term  appears  in  the  v-numerator.  This  is  easily  seen  when 
one  applies  Cramer's  rule  to  form  the  transfer  functions  9^.  This  same 
observation  can  be  made  about  each  of  the  other  numerators;  i.e.  no 
w-term  in  the  w-numerator. 


(9) 


Naval  Coastal  Systems  Laboratory  Technical  Report  287-76,  Development 
of  the  Equations  of  Motion  and  Transfer  Functions  for  Underwater 
Vehicles,  by  D.  E.  Humphreys,  July  1976. 


TABLE  15 

ROOT  SENSITIVITY  ANALYSES  FOR  DERIVATIVE  VARIATION  OF  MOOX 
(Sheet  1 of  2 ) 


LATERAL  DERIVATIVE  ROOT  VARIATION  PERCENTAGES 


IMAGINARY  ROOT  VARIATIONS:  t -O.ld  to  0.10;  * -0.20  to  0.09;  ft  0.15  to  -0.15} 


**  0.19  to  -0.19;  t+t  -0.003  to  0.003;  A**  0 
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TABLE  15 
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LATEAAL  DERIVATIVE  ROOT  VARIATION  PERCENTAGES 


~r 

CHAR.  EQUATION 

V NUMERATOR 

T NUMERATOR 

♦ NUMERATOR 

y'p 

neg. 

Derivative  Variation 
-100»  +100t 

Derivative  Variation 

-toot  +100* 

Derivative  Variation 

-loot  +100t 

Derivative  Variation 
-lOOt  -MOOt 

t -1.00*  1.00* 

t -I.OOt  I.OOt 

0.03t  -0.03S 

-0.06t  0.06t 

0 0 

* 3- 40*  -3.30* 

* 3.60*  -3. 30* 

0.7tt  -0. 70t 

0 0 

# -o.9U  0.90* 

1 -0.9U  0.90t 

-0.07»  0.07* 

0 0 

0 0 

0 0 

k'p 

neg. 

tt  1.1  Ot  -I.OOt 

tt  l.iot  -I.OOt 

0.002$  -0.002t 

0 -O.OOlt 

0 0 

tt  I.OOt  -I.OOt 

tt  I.OOt  -I.OOt 

-0.006t  0.007t 

0 0 

tt  I.OOt  -I.OOt 

tt  I.OOt  -I.OOt 

0 0 

0 0 

0 0 

0 0 

_£ £j 

ttt  0.05t  -0.05t 

ttt  0.05t  -0.05t 

O.Olt  -o.oit 

0 -O.OOlt 

0 0 

***  0.04t  -O.04t 

***  0.04t  -0.04t 

O.Olt  -o.oit 

0 0 

###  -0.03t  0.03t 

###  -0.03t  0.03t 

0 0 

0 0 

0 0 

0 0 

IMAGINARY  ROOT  VARIATIONS:  t 0.2 A to  -0.24;  * 0.13  to  -0.13;  # 0.26  to  -0.26;  tt  0.52  to 


-0.51;  ttt-0.003  to  0.007;  ***  -0.007  to  0.003;  ###  -0.003  to  0.003 


LATERAL  DERIVATIVE  ROOT  VARIATION  PERCENTAGES 


/ 

CHAR.  EQUATION 

V NUMERATOR 

Y NUMERATOR 

♦ NUMERATOR  ! 

Derivative 

Variation 

Derivative 

Varl atlor 

Derivative 

Variation 

Derivative 

Variation 

-100* 

+1G0t 

-100* 

+100t 

-100* 

+100* 

-toot 

+100* 

Y'f. 

t 

-0.05* 

0.05* 

a 

0.06* 

-0.07* 

0 

0 

-0.27* 

0.26* 

t 

-0 . 05t 

0.05* 

* 

0.06* 

-0.07* 

0 

0 

3.70* 

-3.50* 

-0.19t 

0.19* 

-3.90* 

4.20* 

0 

0 

0 

0 

0.15* 

-0.15* 

0 

0 

min. 

0 

0 

K'r 

+t 

-0.06t 

0.06* 

0.16* 

-0.15* 

0 

0 

0.78* 

-0.75* 

t+ 

-0.06* 

0.06* 

AA 

0.16* 

-0.15* 

0 

0 

6.60* 

-5.70* 

-o.oit 

o.oit 

0.03* 

-0.03* 

0 

0 

0 

0 

-0.001 

0 

0 

0 

mi  n. 

0 

0 

N'f 

ttt 

-l.iot 

0.95* 

AAA 

1.30* 

-0.65* 

0 

0 

15.50* 

-10.30* 

ttt 

-1  .lot 

0.95* 

AAA 

1.30* 

-0.65* 

0 

0 

50.20* 

-21.70* 

88.20* 

-31.60* 

96  20* 

-32.80* 

0 

0 

0 

0 

0.53* 

-0.67* 

0 

0 

naj  . 

0 

0 

IMAGINARY  ROOT  VARIATIONS:  t 0.003  to  0;  * -0.003  to  0.003;  ++  0.007  to  -0.007; 


*«  0.003  to  -0.003;  ttt  0.05  to  -0.007;  ***  0.05  to  -0.07 
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APPENDIX  A 

DERIVATION  OF  HYDRODYNAMIC  COEFFICIENTS  RELATIONSHIPS 
FROM  ENERGY  EQUATIONS 


Beginning  with  the  momentum  equations  given  in  Equations  (1)  and 
(2)  in  the  body  of  this  report,  a derivation  is  given  here  for  the  com- 
plete expressions,  showing  the  equalities  and  relationships  Detween  all 
of  the  acceleration  and  nonlinear  hydrodynamic  coefficients  possible 
for  a vehicle  moving  in  an  infinite,  inviscid,  fixed  fluid  with  zero 
circulation. 

A summary  of  the  identities  and  relationships  derived  is  presented 
on  Page  A-6  with  simplifications  arising  from  vehicle  symmetry  on  Page 
A-8. 


Note:  The  hydrodynamic  coefficients  in  parenthesis  with  a subscript 

are  added  together  to  get  the  value  of  the  parenthesized  hydro- 
dynamic  coefficient,  e.g.. 


X 

= (X  ) + (-X  ) 

rq 

rq  ^ rq 

X 

* (M)  + (-N) 

rq 

X 

= M - N . 

rq 

X-FORCE,  SURGE 

. 6T,. 

6T£ 

6T. 

v d f 

f 

f 

"X  ~ dt  Su 

r <5v 

+ q 

6w 

= Au  + B’w  + C'v  + Lp  + Gr  - G’r  + Hq  + H'q 
6T 

-r  -j—  = -Brv  - A’rw  - C'ru  - Mrq  - Frr  - F'rr  - Hrp  + H'rp 

6Tf 

q = Cqw  + A'qv  + B'qu  + Nqr  + Fqq  - F’qq  + Gqp  + G'qp 

letting  Fx  - F + F' , F2  - F - F' , G^  = G + G! , G2  - G - G' , ^ = H + H' , 
H^  = H - H' 


A-l 
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X - -Au  - B'w  - C'v  - Lp  - G2r  - + Brv  + A'rv  + C'ru  + 

Mrq  + Fjrr  * H2rp  - Cqw  - A'qv  - B'qu  - Nqr  - F2qq  - G^p 


A - -X. 
u 


A'  - X = -X 

rw  qv 


L - -X. 
P 


F = X 
1 rr 


F2  ■ -V, 


B = X 

rv 

C = -X 


qw 


B'  = -X.  = -X 

w qu 


C'  = -X.  - X 

v ru 


M = (X  ) G1  a ~X  Qn  m -X. 

rq^  1 qp  2 r 

N - -(X  ) H1  = -X.  H « X 

qr  , 1 q 2 rp 


Other  resulting  relations:  X ■ M - N . 


Y-FORCE,  SWAY 


-Y  = 


dt  Sv 


6T, 


6T, 


^ 6w  + r 6u 


<5T, 


4-  —i  = Bv  + A'w  + C'u  + Mq  + Fr  + F’r  + Hp  - H'p 
dt  <5v 


<5T 

-p  - — = -Cpw  - A'pv  - B'pu  - Npr  - Fpq  + F'pq  - Gpp  - G'pp 

uW 


6T 

r = Aru  + B' rw  + C'rv  + Lrp  - Grr  - G'rr  + Hrq  + H'rq 
ou 


Y = 

-Bv  - 

A'w  - C' 

u - Mq  ■ 

■ V 

- h2p 

+ Cpw 

+ A'pv 

+ B'pw 

+ 

Npr  ■ 

F2pq 

f Gjpp  - 

Aru  - B 

'rw  - 

C' rv  - Lrp  - 

' G2rr  " 

HjTg 

A =■ 

-Y 

ru 

A' 

= -Y.  =* 
w 

Y 

pv 

M = 

-Y. 

q 

F1 

- -Y. 
r 

F2 

= 

Y 

pq 

B = 

-Y. 

V 

B' 

= -Y 

rw 

* Y 

pu 

N = 

(Ynr) 
Pr  x 

G1 

« Y 

PP 

G2 

= 

-Y 

rr 

C - 

Y 

pw 

C' 

- -Y.  - 
u 

-Y 

rv 

L » 

(-Y  ) 
rp 

2 ^ 

- -Y 

rq 

H2 

- 

-Y. 

P 

Y = N - L. 
rp 


A- 2 


Other  resulting  relations: 
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Z-FORCE,  HEAVE 


d ff.  iTf  + 

dt  6w  q *••  ^ 


f A _i 

6u  ^ 6v 


4-  —•  - Cw  + A’v  + B'u  + Nr  + F'q  + Gp  + G’p 
dt  6w 


- q -g—  « -Aqu  - B’qw  - C’qv  - Lqp  - Gqr  + G’qr  - Hqq  - Hfqq 
6Tf 

p -r — » Bpv  + A'pw  + C'pu  + Mpq  + Fpr  + F’pr  + Hpp  - H'pp 


Z = -Cw  - A’v  - B’u  - Nr  - ?2q  - Gjp  + Aqu  + B'qw  + C'qv 
+ Lqp  + G2qr  + H^q  - Bpv  - A'pw  - C'pu  - Mpq  - Fjpr  - H2pp 


A = 

Z 

qu 

A’ 

= -Z.  = -Z 
v pw 

L * (*  ) 
qp  i 

B = 

-Z 

pv 

B’ 

= Z = -Z. 
qw  u 

M » (-Z  ) 

pq 

C = 

-Z. 

w 

C’ 

= Z - -Z 
qv  pu 

N = -Z. 
r 

F.  - -Z 
1 pr 


G.  - -Z. 
1 P 


H « Z 
1 qq 


Other  resulting  relations:  Z =>  L - M. 

pq 


K-MOMENT,  ROLL 

d 6T 

K=- 


6Tf  <5Tf  6Tf  6Tf  6Tf 

K = -d?^r  + w^r-v^r  + r6r_q6r 


j -ST 

6T 

w -g — “ = Bwv  + A’ww  + C’wu  + Mwq  + FjWr  + H2wp 


- v = Cvw  - A’w  - B’vu  - Nvr  - F2vq  - G^vp 


r = Qrq  + P’rr  + R’rp  + Mrv  + F,,rw  + H^ru 


-q  “T — “ 3 Rqr  - P’qq  - Q’qp  - Nqw  - F^v  - G,qu 


F2  - -Z4 


G„  = -Z 
2 qr 


H.  * -Z 

2 pp 
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A*  - K 


ww 


B - (K'  ) 


-K 

L 

■ 

-K. 

F.  - (K 

) - 

(-K 

w 

u 

i wr  ^ 

B'  - -K 

M 

m 

(K  ) - 

(K  j 

G.  - 

K 

vu 

wq  x 

rv  ^ 

1 

vp 

C'  - K 

N 

m 

(-K  ) - 

(-K  ) 

H,  - 

K 

wu 

^ 2 

qw  2 

1 

ru 

(K  ) 

P 

m 

-K. 

i 

X 

i 

-K 

rw  2 

P 

rr 

qq 

Q 

m 

(K  ) 

Q’  - -K.  - 

-K 

rq  1 

r 

qp 

R 

m 

(-K  ) 

R'  -K.  - 

K 

qr72 

q 

rp 

qv 


w 


C - (-K  ) 
vw  , 


F2  ■ <'V. 


G,  - -K 
2 qu 


= -K.  ® K 
2 v wp 


Other  resulting  relations:  K - B - C 


K » M - N 
vr 


K - M - N 
qw 


Krw  “ F1  + F2  " 2F 


K - -P,  - P,  - -2F 
vq  12 


M-MOMENT,  PITCH 

dt  6q 


K *=  Q - R 
qr 


. <5T,  6T  - 6T.  6T,  6T, 

w a i . r f,  t t 

* ',f  u <$w  w 5u  ^ fir  r 


5p 


6T 

-le  «r--(*-p'f  -»’*-»*- v-  v 

6T 

u ■ Cuw  + A'uv  + B'uu  + Nur  + F2uq  + G^up 
6T 

- w ■ -Awu  - B' ww  - C’wv  - Lwp  - G^wr  - H^wq 


A-4 


UtmmaatuuUal. - ...  .. 
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6Tf 

p * Rpr  + P'pq  4-  Q’pp  + Npw  + F^pv  + 

6T. 


-r 


6p 


- Prp  - Q’rr  - R'rq  - Lru  - G.rw  - H_rv 


A - (-M  ) A'  - M 


wu 


uv 


(~M  ) " (~M ) 


wp 


ru 


F.-  M 
1 pv 


B'  - -M  - M 

ww  uu 


M - -Mo 


G.  “ (M  ) 

1 up' 


(-M  ) 
rw  . 


C - (M  ) C'  - -M  N - (M  } « (M  ) 

uw  „ wv  pw'  0 ur', 


- -M. 
1 u 


-M 


wq 


F,  - -M.  - M P - (-M  ) 

2 w uq  rp  . 


P’  - -M.  - M 

r pq 


G2  - <-Mwr)2  - (Mpu)z  Q - -M. 


Q’  = -M  - M 

rr  pp 


H0  - -M 
2 rv 


R - (M  ) 

pr  2 


R’  - -M.  - -M 

P rq 


Other  resulting  relations:  M ■ C - A 

wv 

M - N - L 
pw 


N-MOMENT,  YAW 


M - N - L 
ru 

M - G.  + G0  - 2G 
pu  1 2 

M - -G.  - G « -2G 
rw  1 2 

M - R - P . 
rP 


5T  6T  6T  6T  ST 

N-  - dTsr  + v—  u s^  + « sT- p sr 


A- 5 


m f*  ■'  i 
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j «ST 


<ST, 


« Avu  + B'vw  + C'w  + Lvp  + G2vr  + H^vq 
<ST. 


-u  -r — = -Buv  - A'uw  - C'uu  - Muq  - F.ur  - H0up 
<5v  1 ^ 


q -g — = Pqp  + Q'qr  + P'qq  + Lqu  + G^qw  + H^qv 


6p 

6T 


-p  = - Qpq  - P'pr  - R’pp  - Mpv  - F2pw  - H^pu 


A » 


(N  ) 
vu  . 


A’  = N 


uw 


L - (N  ) - (N  ) 


VP'  i ' W i 


B - (-N  ) 

uv  , 


B’  = N 


vw 


C’  - N * -N 

w uu 


N - -N. 
w 


- -N 

pw 

P - (N  ) 
qp  i 

P* 

- -N.  - 

q 

= -N.  - N 
u vr 

Q “ (~N) 

pq  2 

Q’ 

■ -N.  ■ 
P 

- (-N  ) - 

up  2 

(N  ) R = -N. 

qv  2 r 

R’ 

- N 

qq 

P* 


qr 


F1 = -1 


-1 


„ . (-V2 - <-v>2  Gl . » 


- 0 


Other  Resulting  relations:  N 


A - B 


■aanv-rm\M*7Ti4a^wwrvigRW5,5.,;K,WSJfc-M 
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P - -K,  - 

(-M  ) - 

(N 

£ 

rp  L 

qp 

Q - -M.  - 

(K  ) - 

(-N 

q 

rq  1 

pq 

R - -N.  - 

(-K  ) - 

(M 

r 

q 2 

pr 

-M.  * -N.  « K = -K  ■ M - -N 
r q rr  qq  pq  pr 


q*  - -K.  - -N.  <*  -K  - -M  - M - N 
r p qp  rr  pp  qr 


R'  . -K.  - -M.  - K - -M  - N - -N 

q P rp  rq  qq  pp 


When  xz-plane  symmetry  exists,  then  the  following  coefficients  are  zero; 
A’  , C'  , P' , R’ , L,  M,  N,  G,  G' . 

Table  A1  gives  a summary  of  force  and  moment  coefficients  for  a body 
with  xz-plane  of  symmetry.  When  xy-plane  symmetry  exists,  then  the 
following  coefficients  are  zero;  A’,  B* , P',  Q' , L,  M,  N,  H,  H' . 

When  the  yz-plane  symmetry  exists,  then  the  following  coefficients  are 
zero;  B: , C* , Q' , R’ , L,  M,  N,  F,  F*. 

When  two  planes  of  symmetry  exist,  all  of  the  terms  that  are  zero  in 
the  single  plane  cases  are  zero. 


A- 8 
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TABLE  A1 
(Sheet  2 of  2) 


Coefficient  zero  because  of  xz  symmetry 

**Coeffieient  not  present  for  arbitrary  body  with  no  planes 
of  symmetry 
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APPENDIX  B 

COMP, ARISON  OF  ANALYTICAL  PREDICTIONS  WITH  EXPERIMENTAL  DATA 


In  this  appendix  hydrodynamic  coefficients  of  an  operational  U.S. 

Navy  submersible,  obtained  at  the  David  Taylor  Naval  Ship  Research  and 
Development  Center  are  compared  to  corresponding  data  calculated  by  the 
methods  set  forth  in  the  body  of  this  report.  Table  B1  lists  those  15 
acceleration  coefficients  generally  considered  to  be  important,  their 
experimental  and  calculated  values,  and  the  percent  difference  between 
the  two.  Percent  difference  in  this  table  was  obtained  as  follows: 

Experimental  - Calculated 
Experimental 

As  can  be  seen  in  Table  Bl,  percent  difference  is  significant  in 
only  four  of  the  15  cases,  that  is,  those  above  12  percent. 

Three  of  those  six  coefficients  (M. , Z.,  and  K.)  have  been  shown  in 

w q p 

sensitivity  analyses  illustrated  in  the  body  of  this  report  to  be  of  rela- 
tively minor  significance. 

Some  questions  exist  concerning  the  validity  of  using  the  PMM  ex- 
perimental technique  to  obtain  roll  derivatives,  and  specifically,  in 
this  case,  the  fourth  problem  coefficient  listed  in  table,  . 

It  should  be  noted  that  not  all  of  the  coefficients  obtained  by 
DTNSRDC  were  measured.  Some  by  necessity  were  calculated.  Those  marked 
with  an  asterisk  in  Table  Bl  were  calculated  by  some  DTNSRDC  method  and 
not  by  NCSL  methods. 
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TABLE  B1 

COMPARISON  OF  CALCULATED  VERSUS  EXPERIMENTALLY-DETERMINED 
COEFFICIENTS  OF  AN  OPERATIONAL  SUBMERSIBLE 


Coef . 

Experiment 

Calculation 

% Error 

% 

-.16800E-03 

-.60300E-04 

6A.0 

- - 101 10E-01 

- . 10100E-01 

0.1 

-.*99C 01-03 

-.54880E-03 

10.0 

Zq 

17100E-03 

-.60300E-04 

65.0 

V 

-.33000E-05 

-.21500E-06 

93.0 

Yv 

-.10880E-01 

- . 10001 E-01 

8.0 

Xu 

-.20900E-03 

-.18500E-03 

11.5 

*v 

-.51000E-0* 

-.28500E-0A 

A4.0 

Nr 

-.A8400E-03 

-.53570E-03 

10.7 

A 

*wq 

- . 101 10E-01 

- . 10100E-01 

0.1 

Xvr 

. 10880E-01 

. 10001 E-01 

8.0 

Y * 
wp 

.101 10E-01 

. 10100E-01 

0.1 

N * 

pq 

-.A9600E-03 

-.54900E-03 

10.7 

z * 

vp 

-.10880E-01 

- . 10001 E-01 

8.1 

Mrp 

.A8070E-03 

. 53550E-03 

11.* 
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APPENDIX  C 

DESCRIPTION  OF  TYPICAL  HYDRODYNAMIC  VEHICLE 
USED  FOR  SENSITIVITY  ANALYSIS 


This  appendix  contains  a physical  description  of  the  underwater 
vehicle  used  as  the  subject  of  this  report.  Table  Cl  is  a listing  of 
pertinent  geometric  facts  describing  the  vehicle,  and  Tables  C2  through 
C6  list  calculated  derivatives,  moments  and  symmetry  terms.  Figure  Cl 
is  a set  of  orthographic  projections  illustrating  the  vehicle. 


NCSL  TR-327-78 


TABLE  Cl 

SUBJECT  UNDERWATER  VEHICLE  GEOMETRY 

BARE  HULL 

Vehicle  Total  Length,  Ft . .2353AE  02 

Hull  Maximum  Cross-Sectional  Area,  Ft2  78824 E 01 

Vehicle  Total  Weight,  LBF 79482E  C4 

Vehicle  IXX  Moment  of  Inertia,  Slug-Ft2 82148E  0A 

Vehicle  IYY  Moment  of  Inertia,  Slug-Ft2 21616E  06 

Vehicle  IZZ  Moment  of  Inertia,  Slug-Ft2 2 1 6 1 6 E 06 

Vehicle  Displacement  Weight,  LBF 79A82E  0A 

Displaced  Fluid  IYY  Moment  of  Inertia,  Slug-Ft2  . . . 2 1 6 1 6E  06 

X-Distance  from  Nose  to  C.G.,  Ft.  10391E  02 

X-Distance  from  Nose  to  C.B.,  Ft 1 039 1 E 02 

Z-Distance  from  C.G.  to  C.B.,  Ft -.10000E  00 

Bare  Hull  Wetted  Area,  Ft2 . 1 842 1 E 03 

Bare  Hull  Lateral  Area,  Ft2  .58070E  02 

Bare  Hull  Volume,  Ft3 1241 1 E 03 

TAIL  STRUCTURE 

Qveral 1 

X-Distance  from  Nose  to  A.C.  of  Tail,  Ft 2 1 298 E 02 

Z-Location  of  A.C.  of  Tail  WRT  C.G -.10000E  00 

Sweep  Angle,  Deg 00000E  00 

Taper  Ratio 10000E  01 

Effective  Aspect  Ratio 2A020E  01 

Exposed  Span,  Root  to  Tip,  of  Fins  (each  of  A),  Ft.  .13200E  01 

Horizontal  Tail 

Total  Exposed  Planform  Area  (2  fins),  Ft2  29010E  01 

Dihedral  Angle,  Deg 0000OE  00 

Incidence  Angle,  Deg 00000E  00 

Vertical  Tail 


Total  Exposed  Planform  Area  (each  of  2 fins),  Ft2  . . 1 45 1 0 E 01 
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TABLE  C2 

S N A M E NON-DIMENSIONAL  LONGITUDINAL 
STABILITY  DERIVATIVES 


XU 

m 

-.21A08E-02 

ZU 

.0 

MU 

m 

.0 

xw 

m 

.0 

zw 

- 

-.3R464E-01 

MW 

m 

.25379E-02 

XQ 

m 

1 5275E-03 

ZQ 

- 

-.15275E-01 

MQ 

m 

-.73062E-02 

XTH 

m 

.0 

ZTH 

- 

.0 

MTH 

m 

-.29769E-03 

XUD 

m 

-.62385E-03 

ZUD 

- 

.0 

MUD 

m 

.0 

XWD 

m 

.0 

ZWD 

- 

18182E-01 

MWD 

m 

-.1  WOE-03 

XQD 

m 

.0 

ZQD 

- 

l43^0E-03 

MQD 

m 

-.83506E-03 

XX 

m 

.0 

ZX 

m 

.0 

MX 

m 

.0 

XZ 

m 

.0 

11 

c M 

.0 

MZ 

m 

.0 

XDELT 

m 

57867E-03 

ZDELT 

« 

-.57867E-02 

MDELT 

m 

-.28238E-02 

M - . 19043E-01 

I Y - 

.93512E-03 

i 

! 

I 
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TABLE  C3 

NONLINEAR  X-Z  SYMMETRY  TERMS 


TABLE  C4 

SNAKE  NON-DIMENSIONAL  LATERAL 
STABILITY  DERIVATIVES 


Y V 

-.33046F-01 

KV 

■ 

-.15537E-03 

NV 

-.22682E-02 

YP 

17294E-03 

KP 

m 

-.13896E-OA 

NP 

.80153E-04 

YR 

. 1 5280E-01 

KR 

•53A35E-OR 

NR 

-.73083E-02 

YPHI 

.0 

KPHI 

3 

-.29769E-03 

NPHI 

.0 

YPSI 

.0 

KPSI 

m 

,0 

NPSI 

.0 

YVD 

18182E-01 

KVD 

m 

-.77258E-0A 

NVD 

.1A350E-03 

YPD 

.77258E-04 

KPO 

m 

-.3721  IE-06 

NPD 

.60977E-06 

YRD 

. 1A350E-03 

KRD 

m 

-.60977E-06 

NRD 

-.8351  IE-03 

YY 

.0 

KY 

m 

.0 

NY 

.0 

YDELT 

. 13285E-01 

KDELT 

m 

.56A51E-OA 

NDELT 

-.64829E-02 

C-4 
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TABLE  C5 

MASS  MOMENTS  OF  INERTIA 


M - .190A3E-01 

X 

ISI 

1 

.0 

IX  - .35537E-04 

IZ  - 

.93512E-03 

TABLE  C6 

NONLINEAR  X-Z  SYMMETRY  TERMS 


YUU  - .0 

KUU 

" .0 

NUU  - .0 

YWP  - .18182E-01 

KWP 

- -.77258E-0A 

NWP  - -.14340E-03 

KQR 

- -.45882E-07 

KVW 

- .21335E-06 

YPQ.  - .1A3A0E-03  NPQ  - -.83469E-03 
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APPENDIX  D . 

DESCRIPTION  OP  MANEUVERS  USED  IN  TRAJECTORY  SIMULATION 


Figures  D1  through  D4  are  time  history  plots  of  the  four  basic 
maneuvers  used  with  the  subject  underwater  vehicle  to  generate  tra- 
jectory simulation  and  percent  contribution  plots  of  the  various  hydro- 
dynamic  coefficients  discussed  in  the  body  of  this  report.  These 
maneuvers  were  performed  at  a simulated  vehicle  speed  of  14.35  feet  per 
second . 

The  following  definitions  ere  provided  to  assist  the  reader  in  the 
interpretation  of  Figures  D1  through  D4. 

DELS  ■ Stern  plane  deflection,  in  degrees. 

PITCH  ■ Pitch  angle,  in  degrees. 

Z - Depth,  in  feet. 

ALPHA  ■ Angle  of  attack,  in  degrees. 

DELR  - Rudder  deflection,  in  degrees. 

BETA  - Sideslip  angle,  in  degrees. 

R - Yaw  rate,  in  degrees  per  second. 

PHI  - Roll  angle,  in  degrees. 

Y ■ Off-track  distance,  in  feet. 

X ■ Along-track  distance,  in  feet. 


D-l 


FIGURE  Dl.  5 *±30°  DIVE/CLIMB  MANEUVER  TIME  HISTORY  PLOT 
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